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ABSTRACT 

We simulate the effect of latitudinal variations in the location of star spots, 
as well as their magnetic field strength, on stellar angular momentum loss to the 
stellar wind. We use the Michigan solar corona global MagnetoHydroDynamic 
model, which incorporates realistic relation between the magnetic field topology 
and the wind distribution. We find that the spots' location significantly affects 
the stellar wind structure, and as a result, the total mass loss rate and angular 
momentum loss rate. In particular, we find that the angular momentum loss 
rate is controlled by the mass flux when spots are located at low latitudes but 
is controlled by an increased plasma density between the stellar surface and the 
Alfven surface when spots are located at high latitudes. Our results suggest that 
there might be a feedback mechanism between the magnetic field distribution, 
wind distribution, angular momentum loss through the wind, and the motions 
at the convection zone that generate the magnetic field. This feedback might 
explain the role of coronal magnetic fields in stellar dynamos. 

Subject headings: MHD: Stellar Coronae 



1. INTRODUCTION 

Angular Momentum Loss (AML) through 
a magnetized wind can play a significant 
role in stellar evolution. It has been used 
to explain a number of stellar phenom- 
ena such as the evolution of close binaries 
and cataclysmic variables (CVs), evolution 
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of stellar rota tion rates, and dynamo sat 
uration (e.g. 



Ivanova fc Taam 2003. and 



references therein). 

The principle mechanism by which the 
wind carries out angular momentum from 
the rotating object is as follows. The mag- 
netic field is tangled to the surface of the 
star at one end but its other end is opened 
by the frozen-in wind at the point where 
the flow becomes super-Alfvenic (here- 
after referred to in three-dimensions as 
the 'Alfven surface'). A braking torque 
is applied because there is corotation out 
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to some effective radius larger than the 
stellar radi us and the n et result is stellar 
spin-down ( )Parkerlll958l ). 

Another factor affecting the AML in 
this process is the effect of the magnetic 
field on the stellar wind. In the solar 
case, the magnetic field determines not 
only the large scale distribution of the 
wind (blowing along open field lines), but 
it also determines its bi-modal structure 



where Q is the rotation rate, M is the 
mass loss rate, and ta is the radius of the 
Alfven surface. Based on WD's model, a 
large number of more detailed and com- 
plex studies of AML in cool and hot 
st well as in protostellar accre- 



(PhilliDS et al. 


1995; 


Lamers & Cassinelli 


1999; 


Suzuki) 


2006; 


Wane & Sheelev 


2006; 


Cranmer et al. 2007). 


We can then infer 



callv (Mestel 


1968: Holzwarth 20051 and 


numerically f 


KeDDens & Goedbloedl 


2000; 


Matt & Pudritzl 2008a 


bl; Ud-Doula et al. 


2009]). In particular, 


Aibeo et al. ( 


2007) 



that the distribution of stellar winds in 
stars with hot coronae like that of the Sun 
will be tightly related to the coronal mag- 
netic field configuration . As a result, the 
mass loss and AML will also be influenced 
by the magnetic field structure and its sur- 
face distribution. 

Doppler Imaging and Zeeman Doppler 
Imaging (ZDI) observations of a large num- 
ber of stars have revealed that stellar 
magnetic active regions, appear as dark 
spots ('spots' hereafter), can appear in 
all latitudes, including the polar regions 



(Stra^smeierJll996,2001 



have shown that AML can be modified by 
a factor of 2-4 when introducing magnetic 
flux at high latitudes. They explain the 
change in AML cLS cL result of a change in 
the mass flux or the size of the Alfven sur- 
face, but they argue that the change in the 
magnetic flux tubes topology should not 
significantly affect magnetic breaking. 

In all of the studies mentioned above, 
the wind distribution has been specified 
by some ad-hoc assumptions, though with 
a strong physical basis, such as Parker's 
wind solution. Such a "pre-defined" wind 
can be appropriate in the case of an ide- 
alized and simplified magnetic field topol- 



19971 ). There is also some observational 



Donati fc Collier CameroW , in which the steady-state wind struc- 

Even the well 



evidence that the polar regions of active 
stars harbor brigh t coronal structures (e.g. 



Chung et al.ll2004f ) . In order to understand 



AML in stars is is then important to study 
the relation between the distribution of 
stellar magnetic fields and their AML. 
Weber & Davlsl fll967T ) (WD hereafter) 



have shown that for an idealized, spheri- 
cally symmetric configuration, the loss of 
angular momentum by the wind can be 
written as: 



ture can be predicted 
known Parker solution, which provides a 
very good first order approximation for 
the solar wind, is hydrodynamic in na- 
ture and assumes an idealized, spherically 
symmetric field configuration, with the 
flow freely moving along open field lines. 
These studies are very useful for paramet- 
ric study of AML in stars, but they are 
similar to the potential field approximation 



(Altschuler fc NewkirkHl969l: iMcIvor et al. 



2003c iDonatil 120081 ). in which strong as- 
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sumptions about the boundary conditions 
have been made. 
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In order to compute AML in stars more 
realistically, one should calculate the true, 
non-spherical, Alfven surface, as well the 
actual wind distribution that results when 
realistic magnetic fields with higher de- 
gree of complexity (i.e., with star spots) 
are introduced. In this work, we present 
an implementation of a global Magneto- 
HydroDynamic (MHD) model, originally 
developed for the solar corona, to the 
corona of a Sun-like star. This model is 
driven by 'realistic' surface magnetograms 
and it captures the detailed physics of the 
corona. Unlike other models, the wind in 
this model is determined by the magnetic 
field configuration in a self-consistent man- 
ner. We exploit this approach tool to 
investigate the effect of the magnetic field 
topology on the wind structure and AML 
in a more consistent and realistic way. In 
particular, we focus on the effect of varia- 
tions in star spot location and magnetic 
field magnitude on the wind structure, 
mass loss, and AML rates. 

The topic of this paper covers overlap- 
ping issues from both solar and stellar as- 
trophysics, in which different terminology 
is sometime used to describe similar pro- 
cesses, readers unfamiliar with numerical 
models for the solar corona should consult 



Cohen et all (120071 . 120081 ) for a more de- 
tailed description of the terminology em- 
ployed here. 

We describe the numerical model in Sec- 
tion [2] and present numerical tests and re- 
sults in Section [31 In Section H] we dis- 
cuss the implications of the results over the 
long-term in the context of stellar evolu- 
tion, as well as for the short-term varia- 
tions in the stellar magnetic field (the stel- 
lar cycle). We summarize our conclusions 



in Section [51 

2. MODEL DESCRIPTION 

In order to simulate the corona of a Sun- 
like star, we use a global MHD model de- 
veloped for the Solar Corona (SC) at the 
University of Michigan. We describe the 
model here briefly, and r efer the reader to 
Cohen et al.l (120071 . 120081 ) for a more com- 
plete description, supporting observational 
data, referencing, and model validation. 
The model is based on the BATS-R-US 



generic MHD code (IPowell et al.lll9991 ) and 
it is part of the Spac e Weather Model ing 
Framework (SWMF) flToth et al.ll2005[ ). 

The model solves the set of MHD equa- 
tions: 
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which describe the conservation of mass, 
momentum, magnetic flux, and total en- 
ergy. The equations are solved in a Carte- 
sian grid with Adaptive Mesh Refinement 
(AMR) capabilities, while the condition 
V ■ B = is maintained using the Eight 



3 



Wave scheme (IPowell et al.l Il999l ). This 



SC model has been widely used to sim- 
ulate the steady state solar corona, as 





Cohen et al. 


2008; 


Manchester et al. 


2008 


)• 



2.1. Wind Powering in the Numer- 
ical Model 

In order to mimic the stellar wind pow- 
ering, we adopt a non-polytropic approach 
where 7 = C p /C v is not constant. Here C p 
and C v are the specific heats of the gas. 7, 
is observed to be close to unity at the base 
of the solar corona, but it a pproaches the 



value 5/3 above 15 — 20R Q (ITotten et al. 



19951 ). We can assume that the energy 
necessary to power the wind corresponds 
to this change in 7. We can also assume 
that this change in 7 is larger for the fast 
wind and smaller for the slow wind. The 
end result is that with a proper choice of 
the three-dimensional distribution of 7, we 
can mimic the volumetric heating of the 
wind and obtain a steady state MHD so- 
lution for the corona without introducing 
any heating functions. 

In the solar case, we specify the dis- 
tribution of 7 in a way that is consistent 
with observations. By then using the em- 
pirical Wang-Sheeley-Arge (WSA) model 
flWang fc She"elevlll990l : lArge & Pizzo1l2000h , 
we also obtain a wind spatial velocity dis- 
tribution that is consistent with the mag- 
netic field distribution. The WSA model 
relates the magnetic flux tube expansion 
rate, f s , to the final wind speed blowing 
along this tube. The higher the expan- 
sion is, the slower the wind that originates 
from the tube and vice versa. By using the 
WSA model as an input, we constrain our 



model to approach the empirical value of 
the wind speed. 

The WSA model has been developed 
to predict the solar wind speed at 1AU 
based on the potential field approximation 
and magnetogram data. The model con- 
verts the surface magnetic field data into a 
spherical distribution of the terminal so- 
lar wind speed by matching its parame- 
ters with solar wind data at 1 AU. We 
assume that the wind in Sun-like stars is 
not significantly different in its underly- 
ing physics from the solar case, in particu- 
lar the relation between the expansion rate 
and the wind speed. In principle, the de- 
pendency of the wind speed on the expan- 
sion rate, f s , could be fit to stellar wind 
observations such as those presented by 
Wood et all fl2005h . 

We calculate the potential field and 
WSA coefficients in order to obtain the 
spherical distribution of the final wind 
speed. We assume conservation of energy 
along a streamline (Bernoulli Integral) so 
we can relate the value of the speed at one 
end of the field line far from the Sun to 
the value of 7 at the coronal base at the 
other end. If we assume that far from the 
Sun the total energy roughly equals the 
kinetic energy we can obtain the following 
equation: 



70 p GM e 







[7o - 1] po 



Rc. 



(3) 







where G is the gravitational constant, and 
M Q and R Q are the solar mass and radius 
respectively, u wsa is the speed obtained by 
the WSA model and po, po, and 70 are 
the surface values of the pressure, den- 
sity and 7 respectively. Since the bound- 
ary conditions for the pressure and den- 
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sity are assumed to be known, we can ob- 
tain the boundary conditions for 7. In the 
next step, we introduce a radial function 
of 7 so its value changes from the sur- 
face value to a spherically uniform value 
of 5/3 at r = 12.5i? Q . Finally, we solve 
the MHD equations self-consistently with 
non-uniform 7 until the wind solution con- 
verges. 

We stress that this wind solution is de- 
termined by the magnetic field configura- 
tion. Therefore, any change in the field 
topology should change the wind structure 
in the same way that the magnetic field 
dictates the general structure of the coro- 
nal solution. 

3. SIMULATION 

In this work, we focus on the effect of 
star spot latitude on the angular momen- 
tum loss through the stellar wind. We 
have examined five test cases with different 
surface distributions of the magnetic field 
and have calculated the three-dimensional 
steady state solution for each case. The 
first case is a 5G dipolar field with no spots 
at all. The other cases are a 5G dipolar 
field superimposed with five spot pairs on 
each hemisphere, where the spot pairs are 
located at <p s = 36, 108, 180, 252, 324 de- 
grees in longitude. In each test case, we 
modify the spots co-latitude to be 9 S = 
10, 30, 60, 75 degrees. For each case we 
use a spot magnetic field of lkG and then 
repeat the simulation using a 200G field. 
The general description of the magnetic 
field for each latitudinal case can be writ- 
ten as: 

6-<fe) 2 + (0-6> s ) 2 ]/2.5 2 cos( 



where Bd is the 5G dipolar field, A is 
the field strength required for a magnetic 
flux of lkG (or 200G) in a spot area of 
6.25 cos 6 degrees, and the summation is 
over the longitudes, S , for a particular 
latitude, 9 S . 9 and <fi are both in units 
of degrees. The magnetic flux here corre- 
sponds to the quantity that would be ob- 
tained from stellar magnetogram observa- 
tions using Zeeman-Doppler imaging. 

The polarity of each spot pair compared 
to the large-scale dipole is matched to the 
observed polarity of the emerging flux on 
the Sun (IBabcockl Il96ll ). This perfect 
alignment of the active region is proba- 
bly not realistic, since differential rotation 
will cause one side of the spot to trail af- 
ter the other one. Therefore, our three- 
dimensional simulation could be thought 
as a series of two-dimensional axial simu- 
lations. 

In each case, the artificial magnetogram 
and potential field are calculated. Table 
[1] and Figure [1] show the (dipole) com- 
ponents for each case. It is clear that the 
dipole moment gets weaker for low latitude 
spots (less than 30 degrees) while it gets 
stronger with high latitude spots. The po- 
tential field distribution serves as the ini- 
tial condition for the magnetic field in our 
simulation. We then calculate the WSA 
speed distribution and the distribution of 
7 for each case. Finally, we run the MHD 
code until a steady-state is obtained. 

In all test cases, initial and boundary 
condition are determined in the manner de- 



B 



B r , 



-5> 



(4) 



scribed by ICohen et al.l (120071 ). We use a 
Cartesian grid with 10 levels of mesh re- 
finement so that the grid cell size near the 
inner boundary is of the order of 1/200 .R*, 
where R* is the stellar radius. This high 
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resolution is required in order to resolve 
the active region with at least several grid 
cells. We also dynamically refine the grid 
around current sheets that appear in the 
simulation. All MHD computations were 
performed using the Pleiades cluster at 
NASA's NAS center. 

3.1. Wind dependency on the Mag- 
netic Field Structure 

Figure [2] shows the surface radial mag- 
netic field strength for each test case (left) 
and the resulting wind speed distribution 
calculated from the WSA model (right). 
It is obvious here that there is a strong 
dependency of the wind structure on the 
magnetic field topology. The main effect 
in the case of introducing strong spots on 
top of a dipolar field is that the fast wind 
gets eliminated as the spots migrate to- 
wards the pole. In the polar spot case, 
there is no fast wind at all except for two 
strips of moderate wind (500 km s _1 ) at 
mid-latit udes, which is consi stent with the 
study of iMdvor etafl fl2004 ). 

The reason for this behavior can be 
found by observing the detailed struc- 
ture of the magnetic field in the vicin- 
ity of a star spot. Figure [3] shows the 
magnetic field topology near the star 
spot for the equatorial case (left) and 
the polar case (right). The surface is 
at r = 1.04R*, color contours are of ra- 
dial field strength and streamlines repre- 
sent the three-dimensional magnetic field 
lines. The fast wind comes from open field 
regions, where the expansion of the flux 
tubes is small (IWang fc Sheeleylll990l ). In 
the equatorial case, the star spots inter- 
act mostly with the closed streamers and 
only slightly with open flux located at the 



coronal hole boundary. Therefore, there 
will be some closing of open field lines and 
a slight increase in the expansion of the 
tube, resulting in some removal of the fast 
wind. In the polar case, the star spots in- 
teract completely with field lines originally 
open, so that there are no longer regions 
with non-expanding flux tubes and the fast 
wind is completely eliminated. The regions 
with the minimum expansion are located 
at mid-latitude where we find moderate 
wind strips. 

3.2. Mass and Angular Momentum 
Loss Rates 

Once a steady state is obtained for each 
case, we calculate the Alfvenic Mach num- 
ber, Ma = u sw /ua, with ua = B / \fAixp 
being the Alfven speed, and compute the 
Alfven surface as that where Ma = 1- The 
mass loss rate through the wind, M, can 
be calculated by integrating the mass flux 
through the Alfven surface: 



M 



pu ■ d&A, 



(5) 



where da^ represents a surface element on 
the Alfven surface. The mass flux does not 
depend of course on the particular surface 
of integration, since it is conserved. 

In a similar manner, the AML rate, J, 
can be obtained using Eq. [TJ 



J 



n sin 9r A pu ■ da a, (6) 



where ta and 9 are the local radius and 
latitude at each element of the Alfven sur- 
face. 

The accuracy of the results depends 
mostly on the grid resolution. Since we 
use a non-uniform grid, we estimate the 
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accuracy of our calculation for the largest 
grid cells on the Alfven surface, assuming 
the result would be more accurate were 
smaller grid cells used. We found that for 
the largest grid cells AM / M < 0.001 and 
AJ/j < 0.1. 

3.3. Results for the Case of lkG 
Spots 

Figure H] shows the y = plane from the 
three-dimensional steady state solution for 
each test case (top to bottom). The color 
contours are of number density (left) and 
radial wind speed (right). Black stream- 
lines in each panel represent the magnetic 
field lines, while the solid white line rep- 
resents the Alfven surface. The left-right 
asymmetry is due to the fact that the cut 
goes through the spots on the left side, 
while it goes between two spots on the 
right side. 

The existence of the strong magnetic 
field from the spots locates the Alfven sur- 
face in this direction further out at larger 
radial distance than in the pure dipolar 
field case. Another consequence of the in- 
crease in magnetic field strength due to the 
strong spots is that the current sheet is 
wider than in the dipolar case. 

The results of the MHD solution are 
consistent with the input speed distribu- 
tion obtained by the WSA model, and 
the figures shows how the fast wind is be- 
ing eliminated as the spots migrate to- 
wards the pole. As the spots move towards 
the pole, the Alfven surface shrinks and 
is smallest for the case of 60 degrees co- 
latitude. The Alfven surface gets bigger for 
the polar case, since the field strength at 
high latitudes on the surface is higher due 
to the polar spots. In addition, the overall 



slowdown in the wind for more poleward 
spots causes the density to increase, so we 
find a significant increase in the total mass 
within the co-rotating volume between the 
surface of the star and the Alfven surface. 

3.4. Results for the Case of 200G 
Spots 

Figure [5] shows the simulation results 
with spot magnetic fields of 200G. It can 
be seen for this weaker field that the effect 
of fast wind elimination is reduced due to 
the correspondingly weaker interaction of 
the spot and ambient fields. Unlike the 
lkG case for 60 and 75 degrees, in which 
wind over 300 — 400 km s -1 has been com- 
pletely eliminated, in the 200G case there 
is still some faster wind at higher latitudes. 

4. DISCUSSION 

4.1. The Effect of Wind Distribu- 
tion on AML 

Figure [6] shows line plots of mass loss 
rate (top) and angular momentum loss rate 
(bottom), denoted by dM and dJ, as a 
function of spot co-latitude for the lkG 
case (squares) and the 200G case (trian- 
gles) . The rates are normalized to the mass 
loss rate dM = 1.5 x 10~ 14 M & yr~ x , and 
the AML rate dJ = 10 30 g cm 2 s~ 2 , of 
the dipolar case. dM is also consistent 
with the solar mass loss rate calculated us- 
ing t ypical solar win d parameters at 1 AU 
fe.g. lGombosil2004! l 

The loss rates depend on the interplay 
between the wind speed, density, and the 
size of the Alfven surface. In principle, 
one should expect the AML rate to be con- 
trolled mainly by the mass flux. However, 
this seems to be true only when the fast 



7 



wind dominates, as in the dipolar case or 
when the spots are near the equator. 

In the dipole field case with no spots, 
the rates depend on the high mass flux 
due to the fast wind at high latitudes, and 
on the relatively high density and stronger 
magnetic field close to the equator. When 
introducing equatorial spots, they interact 
with the closed field streamers so the den- 
sity structure near the current sheet does 
not change much. In this case, the equato- 
rial effects on the loss rates are not signif- 
icantly different and the AML rate is con- 
trolled by the change in the mass flux. 

When the spots are located at higher 
latitudes the amount of fast wind is de- 
creased and the AML is actually affected 
most by the significant increase in den- 
sity in the volume between the stellar sur- 
face and the Alfven surface. It is pos- 
sible that this high density, slow wind 
is related to the accumulation of Hct- 
emitting gas in "slingshot" prominences 



tating stars (Collier Cameron & Robinsoi 




1989; 


Collier Cameron & Woods 1992; Jeffries 


1993: 


Bvrne et al.lll996: Eibe 1998: 


Donati et al. 


2000: 


Jardine & van Balleeooijen 


2005) 





The star tries to co-rotate the plasma that 
exists between its surface and the Alfven 
surface, but the stronger field and higher 
plasma mass density apply more torque 
on the star. The net result is higher loss 
of angular momentum, which in this case 
might be described as AML to the stel- 
lar outer atmosphere rather to the stellar 
wind. The mass flux is higher for the case 
of polar spots with 200G than for just the 
dipole (as seen at the top of Figure [6]), and 
the net AML is still greater for the case of 
lkG spots (as seen at the bottom of Fig- 



ure [6]). Figure [7] shows a cartoon of the 
effect of spot location on AML. 

The overall trend in Figure [5] inverts be- 
low 10 degrees due to the size of the Alfven 
surface. In the case of lkG spots, more 
fast wind is eliminated and the Alfven sur- 
face is expanded, so the total mass flux and 
AML are lower. In the case of 200G spots, 
the Alfven surface is smaller than the case 
with no spots. While in this case the wind 
distribution is similar to the no spots case, 
the Alfven surface is now located in re- 
gions of higher density so the total mass 
flux through the surface is higher. 

4.2. Implications for the Long Term 
Stellar Evolution 

The appearance of polar spots in fast ro- 
tating stars has been explained by the en- 
hancement of Corioli s force acting on the 
emer ging flux tubes (ISchuessler &: Solanki 
19921 ) or alternatively by the behavior of 
the meridional circulation in fast rotat- 



ing stars (jSchrijver fc Title! 120011 ). In both 
cases the rotation rate affects the ma gnetic 
field distribution. lAibeo et al.l (120071 ) have 
argued that re-distribution of coronal flux 
in stars does not significantly change the 
wind efficiency for the magnetic breaking. 
However, their model took into account 
only existing magnetic flux with no addi- 
tion of new emerging flux (i.e. active re- 
gions or spots). Our results clearly show 
that the AML through the wind is affected 
by the location, as well as the strength, 
of the spots. Therefore, we argue that as 
a result of fast-rotating young stars having 
their new magnetic flux appear in high lat- 
itudes, their winds are organized in such a 
way that the AML and the spin down of 
the star are large. As the rotation rate de- 
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creases, the new magnetic flux appears at 
lower latitudes and consequently the wind 
is organized in a way that the AML de- 
creases as well. 

The aim of our study is to begin to 
understand how angular momentum loss 
changes with surface magnetic topology, 
with the ultimate aim of understanding 
the rotational evolution of stars. On 
more magnetically active stars, we ex- 
pect the strength of the spots, as well 
as the strength of the background dipole, 
to be higher than the Sun-like case. The 
stronger surface field will induce a stronger 
stellar wind, so the 'ambient' mass-loss 
and AML rates from these stars would be 
higher, regardless of the magnetic spot dis- 
tribution. Here, we have examined how 
redistributing the surface field affects the 
stellar wind distribution and its AML. 
The dominant effect we found was the re- 
organization of the open/closed field lines 
for the different spot latitudes, and elim- 
ination of the fast wind for high latitude 
spots. As the strength of the dipole is in- 
creased, the influence of the spots on the 
magnetic field distribution will be dimin- 
ished. For a very strong background dipole 
an order of magnitude or more higher than 
considered here, the effect of the spots is 
likely to be quite small. It is not yet known 
how the dipolar field scales with stellar ac- 
tivity in general, and future theoretical in- 
vestigation of the effects of the dipolar field 
strength on AML would be worthwhile. 



Kawalerl (119881 ) developed a general ex- 
pression for the rotational evolution of low- 
mass stars based on the WD formulism 
(Eqn. [I]) and an assumed linear scaling 
of radial and dipolar surface magnetic 
field strength with rotation rate, Bq cx Q. 



This simple type of approach has proven 
reasonably successful at matching stellar 
rotation velocities vs. age for solar-like 



Kawaler 


1988: 


Krishnamurthi et al. 


Sills et al. 


20oq 


). However, the rela- 



tion fails for earlier ages because the pre- 
dicted AML is too high to explain rapid 
rotators with periods of one to several 
days in younger op en cluster stars (e.g . 



Stauffer et al. 19871: Barnes fc Sofia 1996 



Krishnamurthi et al. 1997). To avoid this 



problem iKrishnamurthi et al.l (119971 ) intro- 
duced a "saturation" in the dependence of 
the surface magnetic field on rotation rate: 
beyond some critical period the surface 
field was assumed constant. 

Rotation periods of one to several 
days also represents the regime in which 
high latitude surface spots are observed 
in single and tid ally-locked binaries (e.g. 
Strassmeierll200ll ). It is in this regime that 
the complications of spot latitude on AML 
discussed here will be important. The re- 
sults indicate that a simple saturation of 
the dynamo itself and representation of 
the surface field through a mixture of ra- 
dial and dipole components is likely an 
inadequate description. Such oversimpli- 
fication might also confuse the division 
between wind AML and internal angular 
momentum redistribution that is thought 
significant a t ages younger t han that of 
the Hyades flSills et al.l Eoooh . Disentan- 
gling these effects is likely to require both 
sophisticated dynamo and stellar structure 
models, in addition to coronal MHD wind 
models such as those presented here. How- 
ever, much empirical progress could still 
be made through application of the latter 
models to stellar surface magnetograms for 
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stars with different rotation periods in or- 
der to probe the non-orthogonal interplay 
between magnetic field configuration and 
AML. 

4.3. Implications for the Short Term 
Stellar Evolution and Stellar 
Dynamos 

One can ask whether there might be a 
feedback mechanism between the rotation, 
AML and the magnetic field structure. It 
has been shown that re-distribution of stel- 
lar magnetic fields affects the spin down 
of stars and reduces their necessi t y for 
dynamo saturation ( 
Holzwarth fc Jardine 



1997 



Solanki et al. 
2005^ Our results 



take a step forward in this approach. We 
suggest that variations in the the spot 
latitudesre-distribution of the surface mag- 
netic field can lead to variations in the 
wind distribution (as observed in the solar 
case) and result, the AML is modi- 
fied. If the variations of the AML affect 
the motions generating the new magnetic 
flux (i.e. the convection zone) through 
small variations of the rotation rate, then 
we speculate that this might be the con- 
nection between the corona and the stellar 
cycle. 

Nevertheless, the feedback from the 
coronal magnetic field on the solar dy- 
namo has not been studied in detail yet. 
It is possible to perform angular momen- 
tum and mass loss rates calculations based 
on the same method described in this pa- 
per, but driven by real solar magnetogram 
data. The predicted variations in AML 
and mass-loss sould be compared with 
variations of rotation rates in the convec- 
tion zone measure d by GONG and MDI 
flAntia et al. l l2008l . e.g). 



5. CONCLUSIONS 

We test the effect of star spot latitu- 
dinal location and magnetic field strength 
on stellar AML via the stellar wind. Un- 
like previous studies, we use a model in 
which the wind distribution is dictated by 
the magnetic field structure. This key fea- 
ture is crucial in order to get a more real- 
istic picture of the corona, in particular in 
terms of AML rate. Even though we use a 
solar-like wind distribution in this work, it 
seems likely that the relationship between 
the solar magnetic field expansion and the 
solar wind is universal and can be applied 
to other Sun-like stars. 

We find that the change in magnetic 
field structure significantly affects the wind 
structure and result, the loss of angu- 
lar momentum. In particular, we find that 
the angular momentum loss rate is con- 
trolled by the mass flux when spots are lo- 
cated at low latitudes. When spots are lo- 
cated at high latitudes however, the AML 
rate is controlled by the density increase as 
a result of fast wind elimination. 

This work is a step further in our un- 
derstanding of the role of magnetic fields 
in stellar evolution, stellar dynamos, and 
stellar coronal structure. A further inves- 
tigation of these topics involving MHD so- 
lutions driven by surf ace field maps of real 



stellar system s (e.g . iHussain et al.l 12007 



Gregory et al.l 120081 ) . as well as further 



study of idealized cases such as the work 
presented here, could be of great interest. 
This work also demonstrates how the de- 
tailed knowledge of the physics of the solar 
corona can be implemented for other stars 
as we study the solar-stellar connection. 
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Co-Latitude No Spots 10 30 60 75 

lkG 5.0 3.8 4.4 5.99 7.035 

200(2 5.0 4.76 4.89 5.22 5.4 

Table 1: The g® (dipole) component of the 
magnetic field harmonic expansion for dif- 
ferent spot co-latitudes. 
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Fig. 1. — The g® (dipole) component of the magnetic field harmonic expansion with spot 
strength of 200G (top) amd lkG (bottom) as a function of co-latitude. 
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Fig. 2. — Surface radial field strength with color contours representing B r [G] (left) and 
resulted input speed distribution calculated by the WSA model with color contours repre- 
senting u wsa [km s~ l ] (right). The plots are for (top to bottom) dipolar field with no spots, 
and spots at 10, 30, 60, 75 degrees from the equator. Figures la-lj are available in the online 
version of the Journal. 
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Fig. 3. — Star spots located at 10° from the equator (left) interact mostly with the closed 
field streamers, while polar spots located at 80° (right) interact mostly with open flux. The 
sphere is of r = 1.04R* so B r of the spots is weaker than the surface value. Figures 2a-2b 
are available in the online version of the Journal. 
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Fig. 4. — A y=0 cuts in the three-dimensional steady state solutions contain lkG spots. Top 
to bottom show the no spots case and the cases in which spots are located at co-latitude 
of 10, 30, 60, and 75 degrees. Color contours are of number density (left) and radial wind 
speed (right). Streamlines represent the magnetic field lines, the white solid line represents 
the Alfven surface. Figures 3a-3j are available in the online version of the Journal. 
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Fig. 5. — A y=0 cuts in the three-dimensional steady state solutions contain 200G spots. 
Top to bottom show the cases in which spots are located at co-latitude of 10, 60, and 75 
degrees. Display is the same as Figure [3j figures 4a-4f are available in the online version of 
the Journal. 



6.0 - 
5.0 - 
4.0 - 

o 

5 3.0 - 
E 

2.0 - 
1.0 - 



0.0 



lkG spots 
-A- 200 G spots 




10 20 30 40 50 60 70 80 90 
Spots Co-Latitude 



6.0 n 



5.0 - 



4.0 - 

o 

5 3.0 H 
—> 
■a 

2.0 H 



1.0 - 



0.0 



-■- lkG spots 
-A- 200 G spots 





10 20 30 40 50 60 
Spots Co-Latitude 



70 



80 



— i 

90 



Fig. 6. — Line plot of mass (top) and angular momentum (bottom) loss rates denoted by dM 
and dJ as a function of spots co-latitude for lkG spots (squares) and 200G spots (triangles). 
Rates are normalized to the loss rates of the dipolar case dM and dJo. 
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Fig. 7. — The effect of spot location on the AML. The inner circle represents the star and the 
outer circle represents the Alfven surface with tentative field lines (dashed lines). The left 
panel shows the case of low-latitude spots, in which the AML is controlled by the high mass 
loss through the fast wind at high latitudes. In the right panel, which shows the situation 
for high-latitude strong spots, the fast wind is eliminated and so does the high mass loss at 
high latitudes. As a result, the total density in the volume is high and applies more torque 
on the star due to the imperfect co-rotation and the fact that it is frozen-in to the magnetic 
field. 
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